basic.funk <- function(inoculum, c=10)
{
	library(deSolve)

	fn <- function(t, X, params) {
		with(c(X, params), {
			dT <- b - d*T - e*T*V
			dY <- e*T*V - u*Y
			dV <- p*Y - c*V
			list(c(dT,dY,dV))
		})
	}

	params <- list(
		b=1, d=0.1, e=0.001, u=0.5, p=1000, c=c
	)

	T0 <- params$b / params$d
	X0 <- c(T=T0, Y=0, V=inoculum)

	ode(X0, seq(0,100,0.1), fn, params, method="euler")
}
